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We investigate the problem of finding the optimal convex decomposition of a bipartite quantum state into 
a separable part and a positive remainder, in which the weight of the separable part is maximal. This weight 
is naturally identified with the degree of separability of the state. In a recent work, the problem was solved 
for two-qubit states using semidefinite programming. In this paper, we describe a procedure to obtain the 
optimal decomposition of a bipartite state of any finite dimension via a sequence of semidefinite relaxations. 
The sequence of decompositions thus obtained is shown to converge to the optimal one. This provides, for the 
first time, a systematic method to determine the so-called optimal Lewenstein-Sanpera decomposition of any 
bipartite state. Numerical results are provided to illustrate this procedure, and the special case of rank-2 states 
is also discussed. 
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I. INTRODUCTION 

In recent years, a large amount of effort has been put into 
the study of quantum entanglement, driven in part by the re- 
alization of its enormous potential as a resource in quantum 
information processing QJ. The separability problem has re- 
ceived particular attention — this asks for the determination 
of whether a given state of a composite system is separable 
or not. This basic question remains an open problem, and 
indeed, it has already been shown to be NP-hard [2]. Con- 
sequently, complete operational criteria for separability are 
known only in special cases or low dimensions. This poses 
major problems for the characterization of entanglement. For 
one, it makes the quantification of entanglement an extremely 
difficult task. Many quantities of interest, such as the best sep- 
arable approximation measure [ 3 ], robustness of entanglement 
[4] and geometric measure of entanglement |5|, involve some 
kind of optimization over the set of separable states, which 
cannot be handled easily. However, as the separable states 
form a convex set, the study of entanglement often leads to 
convex optimization problems |6|. One can then benefit from 
the tremendous advances made recently in the field of convex 
optimization, as well as the increasing availability of pow- 
erful computing equipment. In particular, semidefinite pro- 
gramming (SDP) [7 1 has found its way, very naturally, into a 
variety of problems in quantum information theory, and has 
been used in the context of distillable entanglement [8|, com- 
pletely positive maps [9], entanglement witnesses [ 10), unam- 
biguous state discrimination IfTTll . and linear optics quantum 
gates [12|. 

In this paper, we are concerned with the Lewenstein- 
Sanpera decomposition |3] (LSD), which is essentially a con- 
vex decomposition of a mixed bipartite state into a separable 
part and a positive remainder. The optimal LSD has maximal 
weight on the separable part, and it is this particular decompo- 
sition that we are interested in. The matter of finding the opti- 
mal LSD is hardly trivial, and even in the two-qubit case, the 
optimal LSDs were only known for some special states iTPH - 
[T5l . In a recent paper lfl6l . it was shown that the optimal LSD 
problem for two-qubit states can be rephrased as a semidefi- 
nite program, effectively solving the problem for this simplest 
possible composite system. This is possible due to the Peres- 
Horodecki criterion lfT71fT8l . which turns a troublesome sepa- 



rability constraint into a positivity and therefore tractable con- 
straint. Coupled with an analysis of the dual SDP problem, 
optimality conditions characterizing the optimal LSD of two- 
qubit states were derived. For qubit-qutrits, the same analysis 
can be carried out, because the Peres-Horodecki criterion for 
separability still holds. Generalizing this approach, we will 
show that the optimal LSD of any bipartite state of arbitrary fi- 
nite dimension can be accessed via a sequence of semidefinite 
relaxations of the optimization problem. Each step in such a 
sequence involves solving a semidefinite program, for which 
numerous reliable and efficient solvers are available ||T9l423l . 
This procedure utilizes the separability criterion introduced 
by Doherty et al. Il24ll25l . and the complementary one pro- 
vided by Navascues et al. ||26ll27ll . which involve searching 
for symmetric extensions of the state in question. 

The structure of this paper will be as follows: in Section[H] 
we introduce the notion of Lewenstein-Sanpera decomposi- 
tions, and discuss the main features of semidefinite program- 
ming. In Section III we present our main result — a system- 
atic procedure for obtaining the optimal LSD of an arbitrary 
bipartite state, using the tools of semidefinite programming 
and the symmetric extensions criterion for separability. Re- 
sults from our numerical implementation will be used to il- 
lustrate the convergence properties of our scheme. Finally, in 



Section IV we will explore the optimal LSDs of qubit-qubit 
and qubit-qutrit states, and then explain how one can obtain 
the optimal LSD of arbitrary rank-2 states analytically. 



H. REVIEW OF LEWENSTEIN-SANPERA 

DECOMPOSITIONS AND SEMIDEFINITE 

PROGRAMMING 

A. Lewenstein-Sanpera decompositions 

Given an arbitrary bipartite quantum state p e SCHa ® < Hb), 
one can look for convex decompositions of p into a separable 
part and a positive remainder. Such decompositions are called 
Lewenstein-Sanpera decompositions [3]. Since the set of sep- 
arable states is compact, an optimal decomposition in which 
the separable part has maximal weight certainly exists. Fur- 
thermore, it was shown by Karnas and Lewenstein that this 
optimal decomposition is unique 1 28 1 for systems of any fi- 



nite dimension. We denote the optimal LSD of p by 



A. Symmetric extensions criterion for separability 



P - <Sg se? + (1 - S)SeM = Psep + Sent- 



(1) 



In the previous equation and in what follows, calligraphic font 
is used to indicate quantities that are optimal, while a tilde 
above an operator indicates that it is not normalized to unit 
trace. It is natural to identify the maximal weight S as the de- 
gree of separability of the state p. Finding the optimal LSD for 
an arbitrary state is in fact a convex constrained optimization 
problem, in which one maximizes tr{p sep }, a linear objective 
function, over the convex cone of separable linear operators, 
subject to the constraint that the difference p - p sep - <x ent 
remains positive semidefinite. 



B. Semidefinite programming 

The primal semidefinite program [7| has the following 
form: 



subject to F(x) > 0, 



(2) 



where F{x) = Fq + YIJL\ X 'F' anc ' ■* e ^'"- The given inputs for 
the primal problem are (i) the vector c e R m characterizing the 
objective function, and (ii) the m + 1 hermitian n X n matrices 
Fq,Fi,..., F m defining the linear matrix inequality. The dual 
problem associated with Q is 



maximize -\i{FqZ} 
subject to \x\FiZ) — c,, i - 1, . . 
Z>0. 



,m, 



The existence of bound entangled states in higher dimen- 
sions prevents us from directly generalizing the above method. 
Nevertheless, Doherty, Parrilo, and Spedalieri have estab- 
lished a separability criterion (the DPS criterion) that is re- 
lated to positivity [24 25 1. Specifically, a positive linear op- 
erator on < H A ® < H B is separable iff it admits a Bose sym- 
metric extension to any number of copies of 1ri B . A £-Bose 
symmetric extension (fe-BSE) of p is a positive operator p e 
2$CHa ® "ftf*) such that tr B *-i {p} = p, and p is Bose symmet- 
ric, i.e., p(I A ® nj: ymm ) = p, where fl* ymm denotes the projector 
onto the symmetric subspace of < H® k . If in addition, p is PPT 
with respect to some bipartition AB'\B k ~', we call p a PPT- 
BSE of p. For convenience, we will only consider the partition 
AB^ 2 ^\B^ 2 \ and define partial transposition to be transposi- 
tion in the last [k/2] copies of *Hb- A linear operator is also 
separable iff it admits a &-PPT-BSE for all natural numbers k, 
and we shall call this the DPS-PPT criterion. For brevity, S k , 
will be used when making statements about S k and S k con- 
currently, and similarly for expressions such as DPS-(PPT). 

The DPS-(PPT) criterion allows one to describe the set of 
separable states with a countably infinite set of positivity con- 
ditions. Let us define S and S k to be the sets of (unnormal- 
ized) states which admit a fe-BSE and A:-PPT-BSE respectively, 
and let S denote the convex cone of (unnormalized) separable 
states. Then the sets {S*}^ form a nested sequence of convex 



cones, i.e, S ' 2 S 2 2 ... 2 S, and similarly for {S kx °° 



p'k=\ 



This 



follows from the fact that a state with a (k + 1)-(PPT)-BSE 
necessarily has a £-(PPT)-BSE, obtained by tracing over one 



(3) copy of "He- Furthermore, lim^_ 



The dual variable Z = Z' is subject to m equality constraints, 
defined by the F,-s and c,s specified in the primal program, 
in addition to a condition of non-negativity, Z > 0. If there 
is a Z > satisfying the dual constraints and a x such that 
F(x*) > 0, the primal and dual problems are called strictly 
feasible. Under these conditions, the optimal primal and dual 
objective values are equal, the sets of optimal variables are 
non-empty, and furthermore, F(x op t)Z op t = 0. 



m. OPTIMAL LSD VIA SYMMETRIC EXTENSIONS 

The main difficulty in the optimal LSD problem is that one 
does not know how to properly characterize the set of sepa- 
rable states. As already mentioned, for IHa ® 'Kb — C 2 ® C 2 
or C 3 ® C 2 , the separable states are exactly those that remain 
positive under partial transposition (PPT) lfT71[T8ll . so that the 
positivity of p sep and its partial transpose pj p suffice to en- 
sure separability. In these cases, finding the optimal LSD 
amounts precisely to solving a semidefinite program lfT6ll . In 
this section, we will describe how to treat the general case of 
pe&(<H A ®'H B ). 



B. Approximating 5 from the outside 



Now, for each teN, one can ask for the £-(PPT)-optimal 
decomposition p = /LkQk + (1 - <h)$k — Qk + Sk mto a con- 
vex sum of a fc-(PPT)-Bose symmetric extendible state and a 
positive remainder, in which the weight of the former is max- 
imal. This is a convex optimization problem over the set S k , •., 
with positivity constraints, and is in fact a semidefinite pro- 
gram (see Appendix |A] for details). We then have a sequence 
of SDPs (of increasing size), indexed by k, where in each SDP, 
the objective function to be maximized is identical, while the 
feasible sets {S k ,J£lj converge to S. Thus, even though the 
problem of finding the optimal LSD of an arbitrary state is not 
a semidefinite program per se, and the feasible set S is fur- 
thermore difficult to sample directly, the DPS criterion allows 
one to contruct a hierarchy of SDPs that "approximate" the 
actual optimal LSD problem. Specifically, the sequence of k- 
(PPT)-optimal decompositions converges to the true optimal 
LSD. We shall now prove this statement. 



Consider the sequence of ^-optimal decompositions 
P 



Unit trace section 



AlQ\ +(1 ■ 

hQi + (1 ■ 


- h)?2 


hQk + (1 - 


- h)Sk 



(4) 



which, a priori, is not known to converge. Also, each k- 
optimal decomposition might not be unique; we will pick any 
optimal one as a representative. Observe that the sequence 
{/UlJtli is monotonically decreasing, since the AkS are the op- 
timal values of the same objective function on decreasing fea- 
sible sets, and are bounded from below by the true degree of 
separability S. Therefore, At, — > /t M > S by the monotone 
convergence theorem. On the other hand, the sequence {Qk)^ = i 
might not be convergent, but it is contained in the compact 
set of normalized states, which has the Bolzano-Weierstrass 
property. That is, there is at least a subsequence {QkJ™ =l that 
is convergent (to some state g x ). The limiting state of such a 
subsequence must be separable, because the terms in the sub- 
sequence get arbitrarily close to the closed set S . Therefore, 
the subsequence of fc-optimal decompositions 

p - hxQh + (1 ~ foi)$ki 

= ^k 2 Qk 2 + (1 - h 2 )$k 2 



hnQk n +(1 -fa„)Sk, 



— AeoQeo + (1 — /too)<T« 



(5) 



converges to a valid LSD. By the definition of S, /loo < S, but 
we already had S < Aoc earlier on, so that equality must hold. 
Since the optimal LSD is unique, the limiting decomposition 
in Q is in fact the optimal LSD, i.e., q^ - g sep and g^, = <Tent- 
Note that the argument above holds for any convergent subse- 
quence. Now, if a sequence in a compact set has the property 
that every convergent subsequence has the same limit, then 
the sequence itself converges to that same limit. The proof of 
this assertion is by contradiction. Suppose q% -t* g sep . Then 
there exists an e > and a subsequence of states {Qk'„}™ =l such 
that \\Qk> n - £> sep || > e for all n e N. On the other hand, this 
subsequence is still contained in a compact set, and thus has 
its own convergent subsequence which, by hypothesis, should 
converge to £> sep , in contradiction to Qk -+> £> sep . Hence, we 
must have Qk — > g se p- 

Therefore, the sequence of SDPs described above is a 
"good" approximation to the actual optimal LSD problem in 
the sense that the optimal LSD can be provably obtained as 
the limit of the sequence of fc-optimal decompositions. Note 
that the same argument holds for the fc-PPT-optimal decompo- 
sitions. One thus has a systematic numerical scheme to com- 
pute the optimal LSD of an arbitrary bipartite state, with the 




FIG. 1. Geometry of the approximation scheme using the DPS-(PPT) 
criterion. The feasible set in the k-th problem is the intersection of 
the affine cone S k (l) - p with the negative semidefinite cone; among 
these candidates, Qt is the operator with the largest trace (highest 
position in the diagram). Here, p lies inside S! \S 2 , so that the first 
SDP has an optimal objective value A\ = 1, whereas the second SDP 
gives A% < 1 . 



precision limited only by the computational resources avail- 
able. In practice, each Ak provides an upper bound (of increas- 
ing tightness) for S. Ideally, one would like to have lower 
bounds as well, which would be possible if we had a com- 
plementary scheme to approximate S from the inside. Such a 
scheme exists: the so-called DPS* criterion, which was given 
by Navascues, Owari, and Plenio ||26ll27l . 



C. Approximating S from the inside 



Navascues et al. l27l showed that the following convex 



cones, formed by a suitable perturbation of the sets S k 



S* k = 



(/»' 



— — o- + — B —tx B {o-} ® f- : a g S k , 
k + d B k + cLb ub 



S; k = \(1- e k )o- + e k tr B {cr) ® j- : cr e S k p 



(6) 



satisfy S* k p) c S for all k, and limS*^ = S. Here, d B = 
dimCHg), Ig is the identity operator on < K B , and e^ is defined 

as 

e k = ^^ min (l - x : P^f mod2) W = o} , (7) 



2(d B - 1) 



where P„ is a Jacobi polynomial J29l . In other words, the 
sequence of sets {SfMfLj approximates S from the inside, 
with (the closure of) the limiting set precisely equal to the 



set of separable states. Note, however, that we no longer have which means that the £* -optimal weight At must obey S - e < 
a hierarchy here — the sequence {£**■)}£!] is not an increas- A* k < S for all large k. In other words, we have A* k — > S. 

By repeating the earlier argument for {Qk)™ =v one can easily 
show that every convergent subsequence of {^}^L[ converges 
to £> sep , from which it follows that the entire sequence of k*- 
optimal decompositions converges to the optimal LSD. 

The argument above fails if p does not have full rank, be- 
cause then, the term in square brackets in Q is not decidedly 
positive. One could of course restrict the choice of 6p to lie in 
the support of p, but with this restriction, g sep +5p might not be 
found inside S *,K for any k. As an extreme example, consider 
a separable p that has a rank smaller than ds- An inspection 
of |6]l reveals that S "A contains only states with ranks that are 
least ds- As a result, the only permissible A:* -decomposition of 
p is the trivial one, i.e., At - for all k. On the other hand, the 
true degree of separability is S - 1 , so we will fail to obtain 
the optimal LSD by taking the limit of the trivial £* -optimal 
decompositions (see also Fig. [4]). 



(P)'k=l 

ing sequence of sets. Nevertheless, we can still construct a 
sequence of SDPs to approximate the optimal LSD problem, 
where this time, optimization is carried out over the sets S*, k 
One then obtains a sequence of £* -optimal decompositions 

p = A\q\ + (1 - A\)g\ 

= 4^2 + (i - 4>«2 



Xjfil + (1 - A* k )s* k 



(8) 



where each state g k e S* ( k In this case, the sequence {/^}£lj 
is not necessarily monotonically increasing, but each A* k does 
provide a lower bound for S. 

As before, we would like to prove that this sequence of de- 
compositions converges to the optimal LSD. This is in fact 
the case for any full-rank state p. The proof requires a differ- 
ent argument from the one used for the fc-optimal decomposi- 
tions. Now, for a full-rank state p, there is a strictly positive 
minimum eigenvalue p m {„, as well as a degree of separability 
>S > fitaia > (since p - p m \ n y > for all normalized states y, 
and in particular, for g sep ). Then for any e > 0, after adding a 
traceless 5g to g sep , the remainder 



Q 

P~(S- e)(£> sep + 6q) = — ^— (p - Sp S ep) 



5 

e 



S(S-e) s 

p 6g 



— <>ent + 



S{S-e) 
p 6g 



(9) 



is positive if the term in square brackets is positive. For in- 
stance, whenever \\6g\\ op < ep m i n , where ||-|| op is the operator 
norm 



||A|| op = sup{||Av||:|M|<l}, 



(10) 



D. Numerical results 




we have 



40 



FIG. 2. Optimal objective values for a generic full-rank state in 
£(C 2 <S>C 2 ), using both the DPS-(PPT) and DPS*-(PPT) criteria. The 
dotted line indicates the best upper bound for <S, which is exact in 
this case. The best lower bound is provided by the DPS* criterion, 
and is within 1CT 3 of the best upper bound. 



P Sg >p (e/WUfl) 

e e 

> P ~ P-mijAB 

>0. 
Here, I AB is the identity operator on ( H A ® < H B . Since S * k 



(p) 



(11) 



S, 



for all sufficiently large k, there will be some p* k € S* (j , that 
is contained in the ep m \ n -neighbourhood of g sep . For these 
P* k = £>se P + Sg k , the inequality |fe|| op < ep^ a holds, so that 
the left-hand side of (|9]l is positive. Every such p* k therefore 
appears in a £*-suboptimal decomposition 



p = (S- e)p\ + &\ 



^.remainder 



(12) 



For our numerical work, we introduce a minor modifica- 
tion to the definition of S ': for partial transposition, we use 
the partition A\B for consistency with the usual notion of PPT 
in bipartite states. S ' is thus the set of PPT states. We first 
test out the formulation described above using a full-rank two- 
qubit state (Fig. [2]i. Since there are no bound entangled states 
in this case (S l p = S 2 p = . . . = S), the SDPs using the DPS- 
PPT criterion give the exact degree of separability from the 
very first level in the hierarchy. The lower bounds provided 
by the DPS* and DPS*-PPT criteria are seen to converge to 
S as the number of extensions of fig increases, with the best 
lower bound within 0.001 of S. Although the convergence 
with the DPS criterion is much slower than with the DPS -PPT 




DPS 



DPS* 



60 




FIG. 3. Optimal objective values for a generic full-rank state in 
S(C 4 ® C 2 ). The PPT criterion already provides what appears to be 
the true degree of separability. As in the C 2 ® C 2 case, the difference 
between the best upper and lower bounds for <S is very small. 



FIG. 5. Optimal objective values for a generic full-rank state in 
S(C 5 ® C 2 ). It appears that good upper and lower bounds for <S can 
be obtained as long as qubit extensions are being considered. The 
DPS-PPT criterion performs best. 
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DPS-PPT 



DPS 




FIG. 4. Optimal objective values for a bound entangled state in 
S(C 4 ® C 2 ) (left). A t = 1 with the DPS-PPT criterion reflects that 
p is PPT. The second test with A 2 < 1 then reveals the (bound) en- 
tanglement of p. Asp has reduced-rank, the optimal objective values 
obtained from the DPS*-(PPT) criteria are all zero, and fail to con- 
verge to S. After mixing a small amount of the identity operator with 
p (right), non-trivial ^'-optimal decompositions are obtained, and A* k 
begins to climb towards S. 



FIG. 6. Optimal objective values for a generic full-rank state in 
£>(C 3 ® C 3 ). The maximum number of qutrit extensions that can 
be considered is limited by the available memory. A sizable gap be- 
tween the best bounds remains after 12 extensions were used. This 
gap is comparable to the gap for qubit extensions at the same stage. 
More significantly, the A k values obtained from the DPS-PPT crite- 
rion does not immediately flatten out after a few extensions. 



criterion, we do obtain some useful information, namely that 
p lies in 5' 6 \5' 7 , since A(, - l,Aj < 1. 

Next, we consider the smallest system in which the PPT 
separability criterion fails. A generic full-rank ququart-qubit 
state is used in Fig. [3] The additional PPT constraint is com- 
putationally more expensive, but convergence is again much 
faster (with respect to the number of extensions used). The 
PPT criterion with a single copy of fis already attains what 
appears to be the true degree of separability. Also, the gap be- 
tween the best upper bound and best lower bound is less than 
0.005. 

We are also interested in what happens for a bound en- 



tangled state. In Fig. HI the state p used for the plot on the 
left is the bound entangled ququart-qubit state from the one- 
parameter family introduced in 1 30], with b — 0.5. Note that 
the first PPT test gives an optimal objective value of 1, consis- 
tent with the fact thatp is PPT-entangled, lying in S „\Sp. The 
optimal objective values obtained from the DPS* approxima- 
tion fail to converge, because p has reduced-rank (see Section 
IIICi. After mixing a small amount of the identity operator 



with p, one still has a bound entangled state, which is now 
full-rank. Then, the DPS* approximation works, albeit rather 
slowly, as is evident from the plot on the right of Fig. [4] 

The procedure continues to perform well for p e S(C 5 ® 
C 2 ) (Fig.pji, where the gap between the best upper bound and 



best lower bound for S is less than 0.001. When we move 
from qubit extensions to qutrit extensions, the dimension of 
the SDP grows faster with k, with computational resources 
limiting the number of extensions that can be considered. For 
instance, in Fig. [6] there is still a gap of about 0.15 between 
the best upper bound and the best lower bound after using up 
to twelve qutrit extensions. 

We remark that in our calculations, memory issues prevent 
us from computing the A:-PPT-optimal decompositions beyond 
a certain value of k. This is reflected in the early termination 
of the PPT curves in Figs. |2]to|6] and will be discussed in the 
next subsection. 



E. Complexity considerations 



is dropped. When memory is a limiting factor, the DPS* crite- 
rion allows more extensions to be accessed, which may some- 
times lead to better lower bounds than those obtained using 
the DPS*-PPT criterion with fewer extensions (see Figs.[2]to 



IV. OPTIMAL LSD OF LOW-RANK STATES 

In this section, we will investigate the optimal LSDs of 
some low-rank states. The first subsection deals with two- 
qubit states, and is followed by a straightforward generaliza- 
tion to qubit-qutrit states. Rank-2 states in any dimension are 
then considered in the last subsection. 



At first glance, the size of the semidefinite program seems 
to increase at least exponentially with the number of exten- 
sions k used. After all, the dimension of "Ha ® f H B k is dAd k B , 
which is exponential in k, and we are presumably optimizing 
over SCHa ® < H B )- However, it actually suffices to consider 



A. Two-qubit states 



positive linear operators acting on the symmetric (with respect 
to interchange ofH B ) subspace ( HA® < J~^\ ymm - The dimension 
of this subspace is 



dim(^ A ®^ ymra ) = ^ 



da + k- I 



d A (k + 1) . . . (d B + k - 1) 

(d B -l)l 
d A (dB + k-\) d «- 1 

(d B -iy. 



(13) 



which is at most polynomial in k. In particular, when d B - 2, 
the growth in dimension is linear, which greatly facilitates 
the large-fc calculations for qudit-qubit states. The number 
of real variables needed to parameterize a candidate p^ is 
dimCHA ® f H% k s mm ) 2 - The numerical algorithms used for 
SDP problems typically involve solving least-squares prob- 
lems, each requiring 0(m 2 n 2 ) steps, where m is the number 
of variables and n is the size of the matrices involved in the 
problem. The number of iterations required scales no worse 
than 0(n^ 2 ). Therefore, for fixed dA,d B , the complexity of 
the SDPs with the DPS criterion scales polynomially with k. 

If one wishes to consider PPT-symmetric extensions, it suf- 
fices to consider positive operators on a subspace of "Ha ®'H B k 
that is isomorphic to <H A ® ?<g™ M^ E3- Impos- 
ing an additional PPT criterion thus introduces an additional 
matrix block of size n' x n', where 



d B + \kjT\ - 1\ ld B + [k/2] - 1 



r*/2i 



lk/2} 



(14) 



so the complexity with the DPS -PPT criterion remains poly- 
nomial. However, the memory resources needed to handle the 
larger matrices when using the DPS-PPT criterion is gener- 
ally much greater than that required when the PPT constraint 



only the subset 8(<H A ® <H™ ; ) of S(«H A ®Hf), namely, the s 



It was also noted B25I that the entire sequence of tests is 
required to fully characterize S for aU x d B > 6. That is, for 
every k, there exist entangled states lying in S k ) but not in 



k+] 



, , Furthermore, for every k, the volume of the set of such 
states, as quantified by the measure introduced in [31 ], is non- 
zero. Only when d A x d B < 6 does a kth PPT test suffice (the 
first one). In these cases, a more detailed study of the opti- 
mal LSDs is possible. For two qubits, analytical expressions 
for the optimal LSDs were found for some special classes of 
states [13 1, while an algebraic way for handling the remain- 
ing states was described in lfT51 . A unifying approach using 
semidefinite programming was recently reported in lfl6l . The 
following is essentially a summary of that work, presented 
here to illustrate the use of SDP in the optimal LSD problem. 

One of the salient features of semidefinite programming is 
its well-developed duality theory, which one can exploit to 
extract valuable information about the problem at hand. For 
example, one can obtain optimality conditions using the so- 
called complementary slackness condition. This says that un- 
der the conditions of strict primal and dual feasibility, the op- 
timal primal and dual variables have orthogonal ranges. 

In the optimal LSD problem for two-qubits, —S is found 
by minimizing the linear objective function -tr{p sep } subject 
to three positivity constraints, namely, p sep > 0,p se B p > 0, and 
p - Psep ^ 0, which we can write as a single matrix inequality 



Psep 













Psep 







[ o 





p 


-(■ 



Psep_ 



>0. 



(15) 



We parameterize p sep by its sixteen real components with 
respect to some basis for S(C 2 ® C 2 ). For example, we 
can use the basis {<r,- ® T/}? _ , where <xq = to = I2 and 
o-j,Tj,i,j - 1,2, 3 are the Pauli matrices. Writing these ba- 
sis elements as {Ej} 1 ^, and defining F-, — dmg(Ej,E j s ,—Ei) 
and c 1 — (-1, 0, . . . , 0), we obtain the primal SDP in the form 
Q. As for the corresponding dual variable Z, we may also 
consider it to be block-diagonal: Z = diag(Zi,Z2,Z^). When 
p is full-rank, one easily verifies that both the primal and dual 



problems are strictly feasible (choose p sep and Z to be suitable 
multiples of the identity). Consequently, we have the neces- 
sary and sufficient condition for optimality: 



(i) 



(ii) 



(iii) 



£?sep 





1 





fsep 











Sent J 


£?sep 











p Jb 

fsep 











Sent. 



\Z\ o 
Zi o 
Zi 



>o, 










Zi 







Z 3 



>0, tr{FiZ] =c h i = 1,...,16. 



(16) 
The dual constraints can be used to express Z,i in terms of 
Z\ and Zi, and after some algebra, we arrive at the optimality 
conditions 



(Zl + Z 2 B )Sent 
elepZl 



~<Tent> 



0. 



(17) 



These are identical to the equations originally found by 
Wellens and Kus [15], here derived using the duality theory of 
SDR With an appropriate reparameterization (see ifTBTl ). the 
rank-3 case can be handled likewise, giving rise to another 
set of optimality conditions, hereafter termed the generalized 
Wellens-Kus equations. 



B. Qubit-qutrit states 

The optimal LSD for qubit-qutrit states, can be found in 
an almost identical fashion. In fact, one needs only to select 
an appropriate basis, for example, by replacing \Tj} 3 _ with 
I3 along with the generators of SU(3). For full-rank states, 
the optimal decomposition again satisfies a set of generalized 
Wellens-Kus equations, which are identical to those for the 
two-qubit case, except that g sep can now be rank-4 or rank-5 
while ^ent can be rank-1 or rank-2. Indeed, the entangled re- 
mainder in the optimal LSD can have rank up to <1a + cLb — 1 
|28| . In analogy to the two-qubit case lfT31l . we remark 
here that when g sep and g s B p have their highest allowed rank 
(6 and 5 respectively), and when g eDt is a pure state, g ent is 
in fact maximally entangled. This can be seen as follows. 
Since g sep has full-rank, Z\ vanishes by the complementary 
slackness condition. Meanwhile, Z.2 must be rank-1 since its 
range is orthogonal to p se s p , which has rank 5. Now, it can be 
shown that for a pure state y e S(C 2 ® C^) with concurrence 
q ll32H34l (normalized to take values in [0, 1]), the partially 
transposed state y TB has eigenvalues -^-,±\, and 0, where 

p = yl - q 2 . In particular, the eigenstate associated with 
the non-degenerate negative eigenvalue is a maximally entan- 
gled state. Therefore, the first eigenvalue equation in ( fT7| i with 
Zi — states that ^ ent is maximally entangled. 



C. Optimal LSD of rank-2 states in C M <g> C" 

The optimal LSD for rank-2 two-qubit states was worked 
out in [14). Here, we describe a prescription to find the op- 
timal LSD for rank-2 states in C M ® C". No use is made of 
SDP, as the optimization can actually be carried out analyt- 
ically. First, we recall some early definitions and theorems 
introduced in |3] [28) . 

Definition 1. A non-negative parameter A is said to be max- 
imal with respect to a positive (possibly unnormalized) oper- 
ator p and a projector P — \ifr) (t//\ if p — hP is positive, but 
p — (A + e)P is not positive for every e > 0. 

Definition 2. A pair of non-negative parameter (Ai,A2) is 
said to be maximal with respect to a positive operator p and 
a pair of projectors (Pi = \ifn) (ifri\ , P2 = tyziifaW if ( i) 
p—(AiPi +A2-P2) is positive, (ii) Aj is maximal with respect to 
p — A2P2 and Pi while A2 is maximal with respect to p — Ai Pi 
and P2 (in the sense of DefinitionU}, (iii) the sum Ai + A2 is 
maximal. 

Maximality as defined above is characterized by the follow- 
ing lemma, proven in (3) ■ Here, Hip) refers to the range of p, 
while p _1 refers to the pseudo-inverse of p if p has reduced- 
rank. 

Lemma 1. The maximal A with respect to p and \\]i) {\[f\ is 
given by 

(a)0,if\i/,)<tK(p). 

(b)(ifr\ p - l \^-\if\^}e'R(p). 

The maximal (Ai,A2) with respect to p and 
CI«Ai > <*Ai I , \4>i) <(fol) is given by 

(a) (0,0), if \if,i), ^2} £K(p), 

(b) (O.^alp -1 l<A 2 r'), if\ijJi)i'R(p)and\^2)e'R(p), 

(c) (((Ailp-'l'AiM^lp-'l'Aa)), 

if\ifri) , \if, 2 ) e K(p) and (^ilp" 1 |<A 2 > = 0, 

(d) ^(mp- l W2)-\mp-'Wi)imp- x Wi) 

-\mp- l \^)\), 

withD = (^ip- 1 hfcX^lp" 1 l-A 2 > - K^iIp" 1 l<A 2 >l 2 , 

if tyi) Mz) e K(p) and (4ri\p- l \if>i),(^2\p- l \^2) > 
K«Ai I p- 1 life) 1*0. 

(e)({i/n\p- l \iff 1 )- 1 ,0), if |(Ai>,|<A 2 > e *<P) and 
((Ailp- 1 |<Ai) > KiAilp- 1 \ifr 2 )\ > (Hp~ x l<A2>. 

The following theorem then characterizes the best separable 
approximation g sep in the optimal LSD J3l l36l . 

Theorem 1. Let V be the set of all normalized product vectors 
v a in the range of p, indexed by a. Then g sep — 2« A a P a iff 
(i) each K a is maximal with respect to p — Y>a'±a h-a'Pa' and 
P a , and (ii) every pair (A a , Ag) is maximal with respect to 
P ~ Za>?a,(s K'Pc' and (P a , Pp). 
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Note that the set of product vectors does not constitute a 
vector space. Furthermore, the set V is generally not finite 
(and indeed, even uncountable), so that direct computation of 
the optimal LSD using Theorem [T] is very difficult. However, 
we will prove that for entangled rank-2 states, V contains at 
most two terms. Then, if V = 0, the optimal LSD is just 
the trivial decomposition. If V — {\if/)}, only condition (i) 
in Theorem [T] is relevant. If V — {|^i),|^2)h there is only 
one pair (Ai, A2) to consider in condition (ii) of Theorem [T] 
which by Definition [2] already implies condition (i). In any 
case, Lemmafllprovides the optimal LSD as it gives explicitly 
the maximal A ff required in £J sep = 2» Ai-^V It remains to 
prove the assertion that an entangled rank-2 state has at most 
two product vectors (up to normalization) in its range. 

In the spectral decomposition of an entangled p ran k-2 - 
A<i \<f>i) (<f>\\ + P2 \4>i) (02l, at least one eigenvector, say \<p\) is 
entangled. Using the Schmidt basis for \<f>\) and assuming that 
M < N,we can write the range of p ra nk-2 as 



^(Prank-2) = span{|0i) , |0 2 >} = span 







Uk 



1M 



Zl2 

Zkk 
ZMM 



(18) 



V ) \Zmnj , 



where the Schmidt coefficients % > appear (in decreasing 
order) in the ((k - l)N + k)-th row, and z,j e C. Now we look 
for product vectors in the range of p ran k-2, i.e., \\j/) — x\<Pi) + 
102) X e C. To facilitate this search, we require a simple 
lemma. 

Lemma 2. A vector a T — (an, ... ,a\N,ai\, ■■■ ,omn) t £ 
C M ® C^, written in a product basis, is a product vector iff 



aiftki — cin^kj for all i,k = 1 , . . . , M and j, I ■ 
Necessity is easy to prove, using 



1, 



,N. 



'a u ' 

flMN, 




■pM, 


® 


, C N, 


= 


'bid' 

b\c N 

b M c\ 

\b M c Nj 



(19) 



{aij ,akj ) - (0,0) but (au,au) + (0,0). But this exceptional 
case cannot happen, because if, for example, a^ + 0, then 



o a kl 



ol a kjo 



would imply that a, ( 



0, which is 



a contradiction. Therefore, (T20D holds for all i, k, I. Next, we 



20|l when we carry out either 
k' . We can thus write V/ in 



note that vf does not depend on i, k, because the same multi- 
plicative factor must appear in 
of the replacements i — > /' or k 
place of vf. Finally, we have 



( ■ \ 
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a n 




«>7o 


'vA 


fl,W 






yN) 


flfci 




a kj a 


'vA 


O-kN 






,vn) 


. : 1 
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1 



(aijA (vA 



\ a Mj ) \Vn) 



(21) 



showing that a is a product vector indeed. 

By Lemma [2j \if/) = x\<pi) + 102) is a product vector iff its 
components with respect to the Schmidt basis of |0i) satisfy a 
set of equations that are at most quadratic in x- Furthermore, 
since at least two Schmidt coefficients are non-zero, there is 
at least one non-trivial quadratic equation in this set, namely 
(CTi + Zn)0cm + Z22) - znzii- Therefore, there can be at 
most two distinct solutions for x> corresponding to at most 
two product vectors in the range of Prank-2- F° r M — N — 2, 
there is only one quadratic equation to be satisfied. It follows 
that there is at least one product vector in the range of p ran k_2, 
which had already been established in [35 ]. For M -2,N -3, 
there are two linear equations and one quadratic equation in x 
that are non-trivial, and the solution set can be empty. This 
means that there are qubit subspaces of C 2 ® C 3 that do not 
contain any product state at all. A rank-2 state with such a 
range space will therefore have only the trivial LSD. 

Note that one does not know a priori whether p lan k-2 is sep- 
arable or entangled. If it is actually separable, two cases can 
arise: (i) the spectral decomposition already gives a decompo- 
sition of Prank-2 into a convex sum of product states, or (ii) at 
least one eigenvector is entangled, and carrying out the above 
procedure gives p sep = p ran k-2- 



V. CONCLUSION AND OUTLOOK 



Now, suppose atjUki = auakj for all i, j, k, I. Since a + 0, there 
must be some io, jo such that a, j + 0. Then a^aa - auakj = 
implies that 



a " I = v ik ai '<> 



V/, k, I, 



(20) 



for some set of complex numbers {vf}, except possibly when 



In this paper, we have expressed the problem of finding the 
optimal LSD of any bipartite state as a sequence of SDPs. 
This allows the efficient numerical computation of optimal 
LSDs. Indeed, we have described how to use the DPS-(PPT) 
separability criterion to form a sequence of semidefinite re- 
laxations to the optimal LSD problem. The corresponding 
sequence of solutions to these SDPs provides a sequence of 
convex decompositions of p which has been shown to con- 



verge to the true optimal LSD. A complementary sequence of 
decompositions is obtained if one uses the DPS*-(PPT) crite- 
rion instead. In this case, we still have convergence provided 
p has full-rank. We have implemented this scheme numeri- 
cally and have found that for qudit-qubit states, the degree of 
separability can be obtained to a good precision with a rea- 
sonable amount of computational resources. The introduc- 
tion of an additional PPT constraint reduces the number of 
extensions that have to be considered, for a given error mar- 
gin, but increases the computational cost significantly. The 
prescription we have provided illustrates the versatility of the 
DPS-criterion — it provides a test for separability, a means to 
compute certain entanglement-related quantities, and as elab- 
orated in this paper, a way to access the optimal LSD of any 
given state. Additionally, we have cast the optimal LSD prob- 
lem for qubit-qutrit states as a SDP, and have shown that the 
remainder in the optimal LSD is maximally entangled in some 
special cases. Finally, we have provided analytically the opti- 
mal LSD of rank-2 states in any dimension. 
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Appendix A: Finding the A -optimal decomposition using SDP 

The k-th optimization problem (without the PPT constraint) 
takes place over the real vector space of hermitian operators 
on ( H A ® ( Hg k . We will denote an operator (not necessarily nor- 
malized) on this space by pk or &k- A feasible pk - tr^t-i {pk\ 
in a valid ^-decomposition p - pk + <x* must satisfy the fol- 
lowing conditions: 

(i) Pk > 0, 



(ii) p-ti B t-i{p k } >0. 



We would also like to highlight a curious link between 
the optimal LSD and the concurrence of a state. For pure 



iM 



-i.<V 



states in C" ® C\ the concurrence is defined l34l l37l by 

C^)) = ^^(l-trfp^}), where p A = tr B {|^><<fl}. This is 
then extended to mixed states via the convex roof construction 
l38l . with C(p) equal to the minimum average pure state con- 
currence, taken over all ensemble decompositions of p. Since 
the concurrence is a convex function which vanishes on sepa- 
rable states, one has the inequality C(p) < (1 — S)C{g enl ). Re- 
markably, equality holds for many classes of two-qubit states, 
for instance, the rank-2 states [ 14 1 and the full-rank states with 
full-rank g sep ifTSI . In other words, the optimal LSD is often 
also an optimal decomposition in the concurrence sense, even 
though the defining properties of the two decompositions are 
quite different. It is not known whether this is a coincidence 
stemming from properties unique to two-qubit states |39|, or 
part of a more general relationship applicable to higher dimen- 
sions. In any case, (1 - S)C(g ent ) provides an upper bound for 
the concurrence of p, which is, in some sense, the best that 
one can do with a decomposition of p into two parts. In fact, 
the bound holds as long as C is a normalized quantity defined 
through a convex-roof construction which vanishes on product 
(and thus separable) states, not necessarily the concurrence. 
Computing convex roof quantities, unfortunately, involves a 
very difficult optimization procedure, so even C(g ent ) is usu- 
ally not accessible. Nevertheless, the quantity 1-5 alone 
is an entanglement monotone in its own right [28 1, and may 
even serve as a fairly good upper bound for the concurrence 
of p. Heuristically, the process of finding the optimal LSD 
"concentrates" all the entanglement properties of p into <r ent , 
whose concurrence is then expected to be close to maximal. 
This suggests that (1 - S) by itself can serve as a fairly good 
upper bound for the concurrence of p. Furthermore, in those 
instances where (1 - S)C(g ent ) - C(p) holds, we even have a 
decomposition realizing the concurrence of p. 



of < H® which has dimension d 



In an actual numerical implementation, we may parameterize 

pk by its components with respect to some orthogonal basis. 

Recall that we only need to consider the symmetric subspace 

d B + k-l\ . ,. 
, . An obvious 

choice would be to use the basis {a-j ® Tj} i=0 t ... x f--ij=Q ,...,cP—i< 

where cr = Id A /dA,To - Id/d, while {o",}/j and {tMj^ are 
the traceless generators of SU(aU) and SU(<i) respectively. In 
this parameterization, pk = Yiu x ijO~i ® Tj< We will also need 
to compute the reduced matrices t, = tr B t-i{Tj}, which enters 
in the expression pk - Y*i,j x ij°~i®Z- Note that among the 
reduced matrices ctj®t ., only <x ® t () has non- vanishing trace 
(equal to 1). Therefore, in standard form, the corresponding 
SDP is: 



subject to F(x) = 0^ ffip 

+ ti j Xij(°~i ® t j) © (-<Ti ® I/) ^ °> 

(Al) 
where Cq = — 1. If the PPT constraint is desired, one addi- 
tionally computes the matrices <x, ® r. sLt/2J and proceeds in 
a similar fashion. For computing the k* -optimal decomposi- 
tions, we replace the reduced matrices cr, ® t . with suitably 
perturbed ones, in accordance with the definitions of S *,K in 



©• 



Appendix B: Optimal LSD from another angle 

There is another way to look at the optimal LSD problem, 
which may be instructive. Let us define a real-valued function 
f p on the set of normalized states, 



f p (y) = supR : p - Xy > 0}, 



(Bl) 
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which takes on values in the closed interval [0,1]. In partic- 
ular, one has f p (g s&p ) - S. The degree of separability of p is 
then the maximal value of f p \ s restricted to the set of separable 
states S . The argument at this maximum, arg max f p \s(y), is 
then p sep . It is possible to show that f p is a quasiconcave func- 
tion, that is, whenever, j\ and 72 are two normalized states, 
and 9 e [0, 1], then 



f p {8y { + (1 - 6) 72 ) > min (/ p ( 7l ), f p (y 2 )) . 
Let us assume that f p {y\) is smaller. Then 



(B2) 



p - /pin) [on + (i - 0)72] = e[p - / P (ri)n] 

(1-0)[p- £(72)72] 

(l~0)[/ p (y 2 )-/ p (ri)]72 



+ ( 
+ ( 



>o, 



(B3) 



from which it follows that f p (0y\ + (1 - ffjyf) > f P {j\), i.e., 
(JB2J holds. Now, an equivalent definition of quasiconcavity 
is that every a-superlevel set {y : f p (y) > a] is convex. This 



gives a nice geometrical picture of our problem using the level 
curves of f p (Fig.fTli. 



fp ~ 12 




FIG. 7. Level curves of f p for a full-rank state with S = 1/3. The 
optimal decomposition is illustrated by the straight line. The sup- 
porting hyperplane of 5 at g sep also supports the f p = 1/3 superlevel 
set. 



[1] R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, 

Rev. Mod. Phys. 81, 865 (2009). 
[2] L. Gurvits, in Proceedings of the Thirty-fifth ACM Symposium 

on Theory of Computing (ACM, New York, 2003), pp. 10-19. 
[3] M. Lewenstein and A. Sanpera, Phys. Rev. Lett. 80, 2261 

(1998). 
[4] G. Vidal and R. Tarrach, Phys. Rev. A 59, 141 (1999). 
[5] T.-C. Wei and P. M. Goldbart, Phys. Rev. A 68, 042307 (2003). 
[6] S. Boyd and L. Vandenberghe, Convex Optimization (Cam- 
bridge University Press, 2004). 
[7] L. Vandenberghe and S. Boyd, SIAM review 38, 49 (1996). 
[8] E. M. Rains, IEEE Trans. Inf. Theory 47, 2921 (2001). 
[9] K. Audenaert and B. D. Moor, Phys. Rev. A 65, 030302 (2002). 
[10] F. G. S. L. Brandao and R. O. Vianna, Phys. Rev. Lett. 93, 

220503 (2004). 
[11] Y. C. Eldar, M. Stojnic and B. Hassibi, Phys. Rev. A 69, 062318 

(2004). 
[12] J. Eisert, Phys. Rev. Lett. 95, 040502 (2005). 
[13] B.-G. Englert and N. Metwally, J. Mod. Opt. 47, 2221 (2000). 
[14] B.-G. Englert and N. Metwally, in Mathematics of Quantum 
Computation, edited by G. Chen and R. K. Brylinski (Chapman 
and Hall, London, 2002), pp. 2575. 
[15] T Wellens and M. Kus, Phys. Rev. A 64, 052302 (2001). 
[16] G. C. Thiang, P. Raynal, and B.-G. Englert, Phys. Rev. A 80, 

052313(2009). 
[17] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 
[18] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 

223, 1 (1996). 
[19] M. Grant and S. Boyd, cvx: matlab software for disciplined 
convex programming, version 1.21 (2010), http://cvxr.com/cvx 
[20] K.-C. Toh, M. J. Todd, and R. J. Tutuncu, sdpt3 4.0 (beta): a 
matlab software for semidefinite-quadratic-linear programming 
(2009), http://www.math.nus.edu.sg/~mattohkc/sdpt3.html 
[21] J. F Sturm, Optim. Methods and Software 11 625 (1999); se- 
dumi: a matlab toolbox for optimization over symmetric cones, 
http://sedumi.ie.lehigh.edu 



[22] J. Lofberg, yalmip: a toolbox for modeling and op- 
timization in matlab, in Proceedings of the CACSD 
Conference, 2004, Taipei, Taiwan (unpublished), 
http://control.ee.ethz.ch/~joloef/yalmip.php 

[23] K. Fujisawa, Y. Futakata, M. Kojima, S. Matsuyama, S. Naka- 
mura, K. Nakata, and M. Yamashita, sdpa-m (SemiDefinite 
Programming Algorithm in matlab Users Manual — Version 
6.2.0, Research Report B-359, Operations Research Depart- 
ment of Mathematical and Computing Sciences, Tokyo Institute 
of Technology, May 2005. 

[24] A. C. Doherty, P. A. Parrilo, and F. M. Spedalieri, Phys. Rev. 
Lett. 88, 187904 (2002). 

[25] A. C. Doherty, P. A. Parrilo, and F M. Spedalieri, Phys. Rev. A 
69, 022308 (2004). 

[26] M. Navascues, M. Owari, and M. B. Plenio, Phys. Rev. Lett. 
103, 160404 (2009). 

[27] M. Navascues, M. Owari, and M. B. Plenio, Phys. Rev. A 80, 
052306 (2009). 

[28] S. Karnas and M. Lewenstein, J. Phys. A 34, 6919 (2001). 

[29] M. Abramowitz and I. A. Stegun, Handbook of Mathematical 
Functions with Formulas, Graphs, and Mathematical Tables, 
(Dover, New York, 1972), pp. 771-802. 

[30] P. Horodecki, Phys. Lett. A, 232, 333 (1997). 

[31] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, 
Phys. Rev. A 58, 883(1998). 

[32] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 

[33] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[34] P. Rungta, V Buzek, C. M. Caves, M. Hillery, and G. J. Mil- 
burn, Phys. Rev. A 64, 042315 (2001). 

[35] A. Sanpera, R. Tarrach, and G. Vidal, Phys. Rev. A 58, 826 
(1998). 

[36] M. Lewenstein, B. Kraus, J. I. Cirac, and P. Horodecki, Phys. 
Rev. A 62, 052310(2000). 

[37] S. Albeverio and S.-M. Fei, J. Opt. B 3, 223 (2001). 

[38] A. Uhlmann, Phys. Rev. A, 62, 032307 (2000). 



11 



[39] K. G. H. Vollbrecht and R. F. Werner, J. Math. Phys. 41, 6772 
(2000). 



